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Abstract. The formation and rupture of atomic-sized contacts is modelled by means of 
molecular dynamics simulations. Such nano-contacts are realized in scanning tunnelling 
microscope and mechanically controlled break junction experiments. These instruments 
routinely measure the conductance across the nano-sized electrodes as they are brought 
into contact and separated, permitting conductance traces to be recorded that are plots 
of conductance versus the distance between the electrodes. One interesting feature of the 
conductance traces is that for some metals and geometric configurations a jump in the value of 
the conductance is observed right before contact between the electrodes, a phenomenon known 
as jump-to-contact. This paper considers, from a computational point of view, the dynamics 
of contact between two gold nano-electrodes. Repeated indentation of the two surfaces on each 
other is performed in two crystallographic orientations of face-centred cubic gold, namely (001) 
and (111). Ultimately, the intention is to identify the structures at the atomic level at the 
moment of first contact between the surfaces, since the value of the conductance is related to 
the minimum cross-section in the contact region. Conductance values obtained in this way are 
determined using first principles electronic transport calculations, with atomic configurations 
taken from the molecular dynamics simulations serving as input structures. 


1. Introduction 

The microelectronics industry faces a great challenge if it is to continue the trend that has 
characterized its success over the last few decades [?]. The size of the gate in a transistor is now 
approaching a few nm in dimension and new strategies must be devised to achieve smaller sizes. 
Nanotechnology is a promising field and potentially offers a solution to the above-mentioned 
problem. In this connection, one of the most important tools both in the development and in 
the characterization of systems at the nano-scale is the scanning tunnelling microscope (STM) 
[?]. With the aid of an STM or mechanically controllable break-junction (MCBJ), another 
related technique, we can study the electronic transport in few-atom systems [?]. Such studies 
typically require two electrodes of the same metal and a piezo system to join and separate 



the electrodes with subatomic precision. By also connecting a battery across the electrodes, 
the current through the nano-contact can be measured, usually in units of the quantum of 
conductance (Go = 2e^//i) [?]. 

The electronic transport through atomic-sized contacts can be interpreted from a plot of 
the measured conductance versus the displacement between the electrodes. Fig. shows 
experimental data for gold electrodes at 4.2K; the blue curve corresponds to separation of 
the electrodes, and the red curve to the process of bringing them into contact. These curves 
are usually referred to as conductance traces. The insets in the figure help to visualise the 
contact formation and rupture processes. Moreover, the traces in the figure exhibit a jump in 
conductance [?] that is close to 1 Gq in value. In this case, the rupture trace corresponds to a 
transition from a monatomic contact to a complete loss of contact. In the case of the formation 
trace we observe a jump in the conductance that corresponds to a transition from tunnelling 
current to first atomic contact [?]. 



Figure 1. (Colour online) Experimental 
conductance measurements recorded at 4.2 K 
using gold electrodes. The insets illustrate the 
geometry of the electrodes as contact is made 
or broken. 


A disadvantage of the above experimental techniques is the difficulty of imaging the electrodes 
at atomic resolution. Therefore, theoretical modelling can prove useful in obtaining more 
information about the geometry of the electrodes and the mechanism of contact formation. 

2. Theory 

Molecular dynamics (MD) simulations [?] based on empirical potentials can serve to model nano¬ 
contact formation and rupture. However, MD relies on interatomic potentials that are fitted to 
a limited number of parameters of the material. So, doubts remain as to the reliability of the 
potentials when systems are taken to situations far away from the range where the parameters 
have been fitted. Also, there is no description of the electrons in MD and it is thus not possible 
to obtain information such as the conductance of the nano-contacts. 

On the other hand, DFT calculations [?] have the advantage of being able to describe both 
the nuclei and the electrons in the system, with the use of certain approximations, but without 
the need of fitting parameters. Electronic transport calculations based on DFT permit the 
conductance to be obtained to a high level of accuracy in materials such as gold [?]. However, 
these calculations are computationally costly and systems of up to only a few tens of atoms can 
be conveniently handled. 

We use a combination of these two methods to study the dynamical formation and breaking 
of gold nano-contacts and to calculate the conductance through the most probable structures at 
the moment of first contact. 

The MD code LAMMPS [?, ?] is used to perform the dynamical calculations. The embedding 
atom potential [?] for gold developed by Zhou et al [?] serves to model the interactions 
between the atoms. This potential has been fitted to parameters such as lattice constants, 
elastic constants, bulk moduli, vacancy formation energies, sublimation energies, and heats of 
solution [?]. 




Note that it has been necessary to modify the LAMMPS source code to enable continuous 
breaking and making of nano-contacts. The simulation structure is divided along its length 
(z-axis) into evenly-spaced bins via a spatial averaging fix that comes standard with LAMMPS. 
Following Ref. [?], the choice of bin size (height) is the interplanar distance between close-packed 
layers (2.04 A along (001) and 2.3556 A along (111)). Special functions have been added to the 
variable subroutine in LAMMPS to find the bin containing the fewest atoms and the number of 
atoms in it. Each rupture-formation simulation is repeated over twenty cycles and paused every 
10000 fs to determine the least-atom bin number and the number of atoms therein. 

System sizes ranging from 500 up to almost 3000 atoms and oriented along two different 
crystallographic directions, (001) and (111), have been considered. The initial configuration 
consists of a neck type structure (see Fig. [^. In the case of a narrow neck (e.g.. Fig. [^a)), the 
direction of motion during each rupture-formation cycle is reversed when there are a minimum 
of 5 atoms in the least-atom bin. In the case of wide neck (e.g., Fig[^b)), on the other hand, 
the cycles are reversed when the least-atom bin contains 7 atoms. Each rupture or formation 
cycle is continued for a further 1000 fs after these conditions have been met to ensure that 
the contact has actually ruptured or that the least-atom bin indeed contains the predefined 
minimum number of atoms. 



Figure 2. (Colour online) Two of the initial 
configurations used in the simulations, left 
panel a) 525 atoms of gold, right panel b) 
2804 atoms. Both are oriented along the (001) 
crystallographic direction in this case. 


The first two top and bottom layers of all the simulated structures are frozen and displaced 
in opposite directions by 0.0041 A and 0.0047 A every 1000 fs for structures in the (001) and 
(111) crystallographic directions, respectively. The velocity of deformation in simulations (0.41 
and 0.47 m/s along (001) and (111), respectively) are within the limits that can be computed 
by MD - about four orders of magnitude below the speed of sound in gold (3000 m/s) (See 
Ref. [?] and the references therein). In order to keep the temperature constant throughout the 
calculation, a Nose-Hoover thermostat is applied [?, ?]. 

For the formation cycles, a Fortran subroutine is used to calculate the minimum distance 
between the two opposing electrodes. Since the electrodes are allowed to separate to about 5 
A, the atoms belonging to the top fragment are labelled ”1” and those belonging to the bottom 
fragment ”2”. The first two distinctly labelled atoms (”1” and ”2”, respectively) to come within 
a distance halfway between first and second nearest neighbours during contact formation cycles, 
are then judged to have made contact first. 

Subsequently, the conductance of selected structures at the moment of first contact - taken 
from the MD simulation trajectories - are computed by means of the electronic transport code 
ALACANT [?, ?, ?], using a minimal basis set that assigns only 1 orbital to each atom. It is 
well known that conductance values obtained from electronic transport calculations on few-atom 
nano-contacts, are determined by the few atoms in the minimum cross-section of the contact. 
Therefore, in order to save time on DFT calculations, we remove a few of the outer layers from 
each electrode. 





























3. Results and Discussion 

3.1. Molecular Dynamics Simulations 

Figures [^a), (b), (c) and (d) show the MD formation traces of the four structures considered 
here. In all cases we observe that the first few cycles of formation (labelled ”runl”, ”run2”, etc.) 
require more steps before first contact, with the exception of ”runl3” in Fig. I^c). All the traces 
clearly exhibit a jump in distance at the moment of first contact between the electrodes. This 
jump is reminiscent of those observed in experimental conductance traces [?]. We also see that 
after a number of cycles the traces all overlap, indicating that the geometric configuration of 
the contact repeats in a regular fashion. When the traces repeat in this way we can talk about 
a process of mechanical annealing or, in other words, that we now have reproducible contacts 
with sharpened tips. 





Figure 3. (Color online) Traces of distance between electrodes versus time for a system of: (a) 
525 atoms oriented along the (001) direction. Insets a, b and c correspond to the geometries 
of the electrodes at the moment of first contact for formation cycles labelled runl, run2 and 
runlS, respectively; (b) 573 atoms along (111). Insets a, b and c correspond to the first-contact 
geometries of formation cycles runl, run4 and run20, respectively; (c) 2804 atoms along (001). 
Insets a, b and c correspond to the first-contact geometries of formation cycles runl, run3 and 
run20, respectively; (d) 2719 atoms along (111). Insets a, b and c correspond to the first-contact 
geometries of formation cycles run4, run8 and runl2, respectively. 

Traces 1 to 6 in Fig. [^a) (525 atoms, (001) orientation) differ significantly from each other: 
the jump occurs at different time steps, and its height also differs from trace to trace. However, 
after run 7 all traces up to run 20 are nearly the same: the jumps occur at the same time 
step and are nearly equal in height. The insets in Fig. [^a) show the atomic configurations 
right before contact for several traces. In many cases contact is made through a single atom, or 
monomer] in other cases, as in inset c, contact is made through two adjacent atoms, or double 
contact. There is greater diversity among the traces in Fig. [^b) (573 atoms, (111) orientation), 
both in terms of the time taken for the jump to occur and its height. Here, it takes at least 
10 indentation cycles before the traces start to overlap. When the traces are stable, the jump 
to contact occurs through a monomer (inset c). However, there are also cases where the jump 

















































occurs through a dimer (inset a) and more than one atom (inset b). In Fig. |^c) (2804 atoms, 
(001)), the traces overlap after only three cycles. Note that the aspect ratio of this neck (the 
length of the neck divided by the diameter of its minimum cross-section) is lower (^2:1) than 
that of the structure in Fig. [^a) (^5:1). As in Fig. |^b), the traces in Fig. |^d) (2719 atoms, 
(111)) exhibit greater diversity. The most energetically stable surface layers of a face-centred 
cubic (FCC) crystal are oriented along the (111) crystallographic direction, since atoms in the 
exposed layers have more nearest neighbours than in any other direction. Because (111) layers 
form on the oblique faces of the pyramid shaped electrodes that are oriented along (001), e.g. 
in Figs. [^a) and i c), it is easier to ’’cold” anneal surfaces in this direction, requiring less than 
half the number of cycles than is needed to achieve reproducible structures along (111). The 
insets in Fig. |^d) show that contact occurs more often through a vertical dimer than through 
a monomer, and double atomic contacts also occur. 

3.2. Results of the DFT ealeulations on the MD-generated eontaets 

Following Refs. [?, ?], we have selected three representative (111) contact geometries from the 
MD atomic configurations at the moment of first contact: the monomer in inset b and the double 
contact in inset c in Fig. 3(b) and the dimer in inset c in Fig. 3(d). We have calculated the 
conductance of these structures since the conductance of their counterparts along (001) have 
already been calculated elsewhere [?]. The results are shown in Fig. 4 and the conductance 
values are expressed in units of Go- The crystallographic orientation and cycle of formation 
of the three configurations are also indicated, as well as the time step of and distance to first 
contact. 
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Figure 4. (Colour online) Different types of 
contact: monomer, dimer and double contact 
and their corresponding conductances. The 
distance to first contact in each case is also 
shown. 


We should expect the conductance of the monomer to be between that of the dimer and the 
double contact, since the tunnelling current across atoms just above and below the monomer 
ought to be higher than the dimer’s - because of the closer proximity of the atoms around the 
monomer. It appears that neighbouring atoms in the minimum cross-section make a substantial 
contribution to conductance in the form of tunnelling current, which could explain the high 
value of 1.60 Gq for the monomer. The size of the distance to first contact is also consistent 
with the type of contact under consideration. Thus, one would expect the jump to be smallest 
in the case of the dimer - the approach of two single atoms towards each other - since the 
oncoming atoms ’’see” fewer atoms, and neighbours on the same electrode effectively keep them 
from jumping to the opposing electrode for longer than would be expected. This contrasts with 
the case of a single atom feeling the attractive force of more than one atom on the opposing 
tip, e.g., the monomer in Fig. 4. The force of attraction between more than one atom on both 
opposing electrodes, e.g., the double contact in Fig. 4, ought to be the greatest and hence so 
should the jump to contact, as is indeed the case. 










4. Conclusion 

We have shown previously [?, ?] that when repeated indentation is performed on structures 
that are oriented in the (001) direction, traces of minimum cross-section versus time overlap 
perfectly. This has been explained as a process of mechanical annealing of the electrodes that 
results in very sharp and stable tips [?]. Here we observe the same behaviour, repetition of the 
traces after several cycles of indentation, but of the minimum distance between the electrodes vs 
time. Moreover, we show that electrodes oriented along (111) can also be mechanically annealed. 
In contrast to (001) electrodes, the traces of (111) electrodes take longer to repeat in a regular 
way, apparently due to the stability of the exposed (111) surface layers on the oblique faces 
of (001) structures. In the case of (111) electrodes we also find that the final jump to contact 
occurs predominantly through a monomer or a dimer, and to a lesser extent, through a double 
contact. 

The conductance values we obtain are in agreement with previous works [?]. However, in the 
case of the monomer it is higher than expected. Instead of using a minimal one-electron basis 
set, it might be necessary to use an eleven-electron basis set to improve this result [?]. 
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